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Abstract 

Background: Pseudomonas aeruginosa {PA) has a non-clonal, epidemic population with a few widely distributed and 
frequently encountered sequence types (STs) called 'high-risk clusters'. Clinical P. aeruginosa {cWnPA) has been studied in all 
inhabited continents excepted in Africa, where a very few isolates have been analyzed. Here, we characterized a collection 
of cWnPA isolates from four countries of West and Central Africa. 

Methodology: 184 non-redundant isolates of clinP/A from hospitals of Senegal, Ivory Coast, Nigeria, and Central African 
Republic were genotyped by IVILST. We assessed their resistance level to antibiotics by agar diffusion and identified 
the extended-spectrum p-lactamases (ESBLs) and metallo-p-lactamases (MBLs) by sequencing. The population structure of 
the species was determined by a nucleotide-based analysis of the entire PA MLST database and further localized on the 
phylogenetic tree (/) the sequence types (STs) of the present collection, {//) the STs by continents, (iii) ESBL- and MBL- 
producing STs from the IVILST database. 

Principal Findings: We found 80 distinct STs, of which 24 had no relationship with any known STs. 'High-risk' international 
clonal complexes (CC155, CC244, CC235) were frequently found in West and Central Africa. The five VIM-2-producing 
isolates belonged to CC233 and CC244. GES-1 and GES-9 enzymes were produced by one CC235 and one ST1469 isolate, 
respectively. We showed the spread of 'high-risk' international clonal complexes, often described as multidrug-resistant on 
other continents, with a fully susceptible phenotype. The MBL- and ESBL-producing STs were scattered throughout the 
phylogenetic tree and our data suggest a poor association between a continent and a specific phylogroup. 

Conclusions: ESBL- and MBL-encoding genes are borne by both successful international clonal complexes and distinct local 
STs in cWnPA of West and Central Africa. Furthermore, our data suggest that the spread of a ST could be either due to its 
antibiotic resistance or to features independent from the resistance to antibiotics. 
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Introduction 

Pseudomonas aeruginosa is an important opportunistic human 
pathogen causing infection in patients with impaired immune 
systems [1]. Its intrinsic resistance to many classes of antibiotics 
and its capacity to acquire resistance to almost all effective 
antibiotics during treatment render infections witli tliis microor- 
ganism very difficult to treat [2,3]. Resistance to ^-lactams is of 
particular concern in clinical practice. Hence, high-level resistance 
to these compounds is achieved by AmpC cephalosporinase 
overproduction or by the production of acquired (^-lactamases 



with an extended spectrum (i.e., extended-spectrum- fi-lactamases 
[ESBLs], metaUo-|3-lactamases [MBLs], and extended-spectrum 
oxaciUinases) [4]. 

From a limited number of studies often focused on antibiotic- 
resistant isolates, it appears that P. aeruginosa has a non-clonal, 
epidemic population with a few widely distributed and frequently 
encountered sequence types (STs) called 'high-risk clusters' [5-10]. 
The population of clinical P. aeruginosa has been studied in 
Europe, Asia, Oceania, North and South Americas. Although 
West and Central Africa comprises large city agglomerations, little 
is known about the resistance level to antibiotics and the 
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epidemiology of clinical P. aeruginosa. As proof, only 0.5% (8 
isolates out of 1542) of the P. aeruginosa multilocus sequence 
typing (MLST) database concerns African isolates [11]. 

Large-scale epidemiological studies rely on the analysis of 
MLST developed by Curran et al. in 2004. This method is a 
reproducible manner to compare isolates that are temporally and 
geographically distant and handled by different laboratories, 
allowing a better knowledge of the global epidemiology in P. 
aeruginosa [7]. Quantilying the diversity of specific house-keeping 
genes by MLST analyses is a powerful approach for understanding 
the evolution of the core genome and the processes that shape the 
species biodiversity. 

The aim of this study was to characterize a collection of clinical 
P. aeruginosa isolates from four countries of West and Central 
Africa. We assessed the resistance level to antipseudomonal agents 
and identified the P-lactamases with an extended-spectrum in the 
resistant isolates. We also identified all their STs and compared 
their distribution in a global phylogenetic tree built with all the 
currently defined STs described worldwide. 

Material and Methods 

Isolate collection 

P. aeruginosa non-redundant clinical isolates were collected 
between 2002 and 2012 from patients hospitahzed in the 
University Hospital of Fann (Dakar, Senegal, n — 78), in the 
Centre National de Reference des Antibiotiques de Cote d'lvoire 
(Institut Pasteur, Abidjan, Ivory Coast, n — 48), in the Department 
of Medical Microbiology and Parasitology of the University of 
Lagos (Lagos, Nigeria, n = 30), and in the Institut Pasteur of 
Bangui (Bangui, Central African Republic, n = 28) for a total of 
184 isolates. All the isolates were stored in Brain-Heart Infusion - 
Glycerol broth at — 80°C and managed at the Centre de 
Ressources Biologiques Ferdinand Cabanne at the University 
Hospital of Besanijon. Data were entirely anonymised throughout 
the entire study and thus ethical approval was not needed. 
Figure 1 represents the geographical origin of tlu- isolates. All 
isolates were identified by matrix-assisted laser desorption 
ionization-time of flight mass spectrometry (MALDI-TOF MS) 
with a Microflex LT (Bruker Daltonik GmbH, Bremen, Ger- 
many), according to the manufacturer procedures [12]. 

Antimicrobial susceptibility testing and p-lactamase 
identification 

We assessed the activity of 9 antibiotics from four different 
classes (non-carbapenem P-lactams: cefepime, piperacUlin-tazo- 
bactam, ticarcillin, ceftazidime, aztreonam; carbapenems: imipen- 
em; aminoglycosides: tobramycin, amikacin; fluoroquinolones: 
ciprofloxacin) against the P. aeruginosa isolates by the disk 
diffusion method recommended by European Committee on 
Antimicrobial Susceptibflity Testing (EUCAST, 2013 [13]). Here, 
all non-susceptible isolates were considered as resistant. We also 
identified extended-spectrum P-lactamases (ESBLs) and metallo-P- 
lactamases (MBLs) in isolates resistant to ceftazidime by the 
phenotypic method described elsewhere [14]. For isolates consid- 
ered positive by this approach, the involved enzymes were 
identified by PCR and sequencing with primers targeting ESBL- 
and MBL-encoding genes [15]. 

Genotyping by pulsed-field gel electrophoresis (PFGE) 

The clonality of all the isolates was investigated by PFGE, with 
Dral digestion, as previously described [16]. GelCompar software 
was used for cluster analysis (Apphed Maths, Kortrijk, Belgium). 
The Dice correlation coefficients were grouped and the UPGMA 



clustering algorithm was used to depict the groups as a 
dendrogram. According to international recommendations, we 
clustered isolates that give PFGE patterns with £3 fragment 
differences in the same pulsotype (PT) [17]. We further genotyped 
by multi-locus sequence typing one isolate per PT. 

Genotyping by multi-locus sequence typing (MLST) 

MLST was performed according to the protocol of Curran et al. 
modified by van Mansfeld et al. [18]. Nucleotide sequences were 
determined for internal fragments of the acsA, aroE, guaA, mutL, 
nuoD, ppsA, and trpE genes and were compared with sequences in 
the P. aeruginosa MLST website (http://pubmlst.org/ 
paeruginosa) for the assignment of allele numbers and sequence 
types (ST) [11]. Clonal complexes (CCs) were defined with the 
software START2 as a group of STs sharing at least 5 loci [19]. 

Analysis of MLST data 

In order to build a dendrogram with th(^ 1595 STs available at 
the time of the study (including the new ST described in this 
collection), we concatenated the sequences of 7 MLST genes to 
form a 2,882-bp sequences alignment, defining 664 polymorphic 
positions. The best-fit nucleotide substitution model for this data 
was GTR-I-G-I-I, as determined with jModelTest 0.1.1 [20]. We 
used the Pseudomonas fluorescens PfD-1 as the outgroup [21]. 
Maximum likelihood tree was constructed with RAxML 7.2.8 
[22] and visualized with Dendroscope 3.2.10 [23]. In every case, 
1000 bootstrap repetitions gave values above 900 for most 
branches. 

Results 

Genotyping of the collection 

We genotyped all the 184 clinical isolates from the four West 
and Central Africa countries by PFGE. We obtained 107 
different PTs. As PFGE typing is more discriminatory than 
MLST, we assum(;d here that all the isolates that give PTs with ' - 
3 fragment differences shared the same ST (ChoUey et al., 
submitted manuscript). We therefore determined the MLST 
profiles of one isolate for each of the 107 PTs. The ST of all the 
isolates is given in the Table SI. We found 80 distinct STs. 
Among them, 47 were already in the MLST database wliik- the 
other 33 (from ST1437 to ST1493) were new [1 1]. Nine of these 
new STs belonged to widespread clonal complexes: ST 1461 and 
ST1477 belonged to the CC195, the ST1467, ST1479, ST1481, 
ST1482 and ST1484 belonged to the CC244, die ST1466 
belonged to the CC274 and ST1472 belonged to the CC233. 
Besides, 24 new STs were singletons, meaning that they share less 
than 5 loci with known STs. In our collection, we sought to 
identify isolates of successful international STs or CCs (i.e. that 
are found worldwide). We identified from the MLST database 15 
intercontinental STs or CCs (isolated in patients from a3 
continents among Europe, North and South Americas, Asia and 
Oceania): ST17, ST27, CC155, ST179, CC195, CC233, CC235, 
CC244, ST254, CC274, ST277, CC395, ST446, ST560, and 
CC654 [1 1]. These STs and CCs were not specific to the human 
host and have also been described in the environment. Nine of 
these widespread clinical STs or CCs were also present in West 
and Central Africa. Hence, CC155, CC244, CC235, CC195, 
CC274, CC233, CC395, CC654, and ST277 were represented 
by 23, 19, 8, 4, 4, 3, 3, 3, and 1 isolate, respectively in the present 
collection. 
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Dakar (Senegal, n=78) • 

GES-9 in ST1469 (n=l) 



Abidjan (Ivory Coast, n=48) 

VIIVI-2 in CC233 (n=2] 
VIM-2 In CC235 (n=2) 




Lagos (Nigeria, n=30) 
GES-1 in CC235 (n=l) 
y\M-2 in CC233 (n=l) 



Bangui (Central African Republic, n=2S) 



Figure 1 . Map of Africa showing the four cities included in the study, the number of Pseudomonas aeruginosa clinical isolates, and 
the localization of ESBL- and MBL-producing isolates. 

doi:! 0.1 371/journal.pone.01 07008.g001 



Resistance to antibiotics and production of extended- 
spectrum p-lactamases (ESBLs) and metallo-p-lactamases 

(MBLs) 

Table 1 details the resistance rates to antibiotics of the isolates 
broken down by the country of origin. The antibiotic resistance 
profile of all the isolates is detailed in the Table S 1 . The resistance 
rates to the tested antibiotics remained low (<10%) for isolates 
from Senegal and Central African Republic, with the exception of 
the higher resistance rate to ticarcUlin in the latter country (25%). 
This contrasts with the much higher resistance rate of isolates from 
Ivory Coast and Nigeria. In these countries, 23-40% of the isolates 
were resistant to penicillins (ticarciUin and piperaciUin/tazobac- 
tam), 13-30% were resistant to cephalosporins (cefepime and 
ceftazidime), 10-20% were resistant to aminoglycosides (tobra- 
mycin and amikacin), and 27% were resistant to ciprofloxacin. 

Eleven out of the 184 isolates collected were resistant to 
imipenem (Ivory Coast, n = 8; Nigeria, n = 3). We identified the 
metaUo-P-lactamase VIM-2 in 5 of them. The remaining isolates 
did not produce any carbapenemase and were most probably 
resistant to imipenem after mutations in the gene encoding OprD 
or the down-regulation of the production of this porin [24] . Three 
VIM-2-producing isolates belonged to the CC233 (Ivor\' Coast, 
w = 2; Nigeria, n- I) and two belonged to the CC244 (Ivory 
Coast, w = 2). 

Sixteen isolates of the whole collection were resistant to 
ceftazidime. They were isolated mostly in Ivory Coast (w = 7) 
and Nigeria {n = 6). One CC235 isolate from Nigeria produced 
GES-1 and one ST 1469 isolate from Senegal produced GES-9. Of 
note, none of the 28 isolates collected in Central African Republic 
produced |3-lactamase with an extended-spectrum. 

Interestingly, most of the isolates (17 out of 25) that displayed a 
multi-drug resistant phenot)'pe according to the consensus 
definition of Magiorakos et al. (i.e., resistance to S: 1 agent in a 
3 antimicrobial categories) belonged to the intercontinental clonal 



complexes CC233 (n = 3), CC235 (n = 5), CC244 [n = 3), CC393 
(w = 3), and CC654 (n = 3) [25] . These multi-drug resistant isolates 
mostly came from Nigeria (n = 1 1) and Ivory Coast (m = 10). 

Spread of wild-type and antibiotic-resistant STs and CCs 
in West and Central Africa 

Among the collection, seven STs were represented by more 
than 5 isolates. CC155 spread in Senegal with 23 isolates with a 
wild-type susceptibihty phenotype. We retrieved the CC244 in the 
four countries studied. Of the 19 CC244 isolates, 4 were multi- 
drug resistant (among which 2 isolates produced VIM-2) while the 
remaining 15 displayed a wild-type susceptibility phenotype. The 
CC235 was represented by 8 isolates found in all countries except 
Senegal. Most of the CC235 isolates (6 out of 8) were resistant to 
& 1 P-lactam, to a 1 aminoglycoside, and to ciprofloxacin, and one 
Nigerian isolate produced the class A ESBL GES-1. All the ST699 
isolates found in the collection of Nigeria, Senegal, and Ivory 
Coast {n = 1) were susceptible to all the antibiotics tested. The 
seven ST856 isolates were retrieved in Nigeria. We also found 7 
isolates of CC393 coming from Nigeria (w = 6) or Central African 
Republic (« = 1). FinaUy, the 5 ST1464 isolates found in Senegal 
were fuUy susceptible to all the tested antibiotics. 

Population structure of clinical Pseudomonas aeruginosa 
in West and Central Africa 

We then spotted the West and Central African P. aeruginosa 
isolates on the maximum likelihood (ML) tree built from the entire 
P. aeruginosa PubMLST isolate database (1595 STs described on 
July 2013) [11]. The ML tree displayed a bush-like structure 
(Figure 2). Clonal complex is defined as a group of STs in a 
population that share S5 alleles. In the vast majority of the cases, 
the STs grouped in one CC clustered on the ML tree (Figure 3). 
For instance, the STs that belong to the CC235 (e.g., ST235, 
ST533, ST534, ST622, ST227, and ST745) clustered on the ML 
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Number of isolates 
■ 1 isolate 

• 2-5 isolates 

• 6-10 isolates 
S 23 isolates 



^ STs with no ESBL or MBL producers 

^ STs with > 1 isolate producing ESBL or IVIBL 

^ STs of reference strains of P. aeruginosa 



P. fluorescens PfO-1 
(distance, 24.0) 



ST1458 



rST777 

A ST1247 
STH70 ~ 
STM75 ^ ^ ST57 
ST146-LESB5B ST277 
ST3B1^T152 ST170 ©STaeO 



/ 



. ST1469 

O ^ S1 102 4-B1 36-33 

^ ST746 



/ // 



ST253-PAH STZ35-HCG2.5i ST316 
STUW ST1M5 



ST672 

* S112:t9-M18 ^ 
^. ST395 



STe74 ■ 

ST234 STH60 ST514 ST430 ^^j^ ST91G 
ST569 ST1437 ST447 




ST1463 
ST1474 

ST1296 • ^iios * 

ST1468ST1S8RP73* ««Tii7i 511145 ^'^^"^ ST242 
SI14B3 s„„« 9 *iT1471 ., ^^^^^ 

ST282 ^, ST1481 

ST1482 * STS99 ST1473 
ST1479 „ 
1 
ST217 



^ ST213 
^* ST499 
ST1490 



ST1480 



ST5G8 
^ SI38fi-DK2 




Figure 2. Distribution of the STs and CCs of the 184 clinical isolates of Pseudomonas aeruginosa isolated in West and Central Africa 
on a dendrogram built with the data all known STs (n= 1595). STs without any isolates producing ESBL or MBL were represented with red 
spots. STs represented by >1 isolate producing ESBL or MBL were represented with green spots. STs of reference isolates are represented with blue 
diamonds. The dendrogram is based on the alignment of the concatenated sequences of the acsA, aroE, guaA, mutL, nuoD, ppsA and trpE genes 
(forming an artificial 2,882-bp sequence) of the 1595 STs of the P. aeruginosa MLST database in July 2013 [1 1]. See Material and methods section for 
details. 

doi:1 0.1 371 /journal.pone.01 07008.g002 



tree. In two cases, STs from the same CC were scattered 
throughout the ML tree. Hence, ST743 and STS 1 1 were distant 
from their respective predicted founders ST233 and ST 155 
(Figure 3). 

The 80 STs and CCs found in the four African countries (red 
dots on tlie Figure 2) were scattered throughout the tree with no 
particular association with a specific phylogroup. The STs or CCs 
that included ESBL- or MBL-producers (CC233, CC235, CC244, 
and ST 1469) also were scattered throughout the ML tree (green 
dots on the Figure 2). We also repeated this experiment with the 
clinical STs and CCs isolated from the 5 other inhabited 
continents (Figure 4), Here again, the STs and CCs from a given 
continent were scattered throughout the global ML tree. 

Discussion 

Here we analyzed the genetic diversity and the population 
structure of 184 clinical P. aeruginosa isolates found in West and 
Central Africa and identified the STs that produced P-lactamases 
with an extended-spectrum (ESBL and MBL). We found that P- 
lactamases with an extended-spectrum in this region are mosdy 
borne by successful international clonal complexes but also by 
distinct local STs. Our data also suggest a very limited association 
between a genotype and a region. 

Epidemiology of ESBL- and MBL-producers 

VIM-2 was the only carbapenemase found in our collection, 
confirming that this enzyme is the dominant metaUo-P-lactamase 
over the world and that it is well estabhshed in the six populated 
continents [26]. Here, the 6ZaviM-2 was borne by the international 



complexes CC233 (n = 3) and CC244 (n = 2). The presence of 
VIM-2 has been established in Africa. Hence, a ST 1488 isolate 
harboring 6/flviM-2 has been reported in Abidjan (Ivory Coast) 
[27]. ST1488 belongs to the CC308, already associated with 
blttMBL [5] ■ Besides, a CC233 isolate producing VIM-2 has already 
been described in Norway and associated with repatriation from 
Ghana [28]. In our collection, the three CC233 isolates (obtained 
from Nigeria and Ivory Coast) produced VIM-2. This confirms 
the presence of the VIM-2 enzyme in western African countries of 
the Gulf of Guinea. The intercontinental CC244 that has been 
described in Europe, Russia, China, Australia, and Brazil has only 
been associated with 6/apER-i in Poland [1 1,29]. It is worth noting 
that we found CC244 isolates in the four African countries and 
that the majority of these isolates (15 out of 19) were fully 
susceptible to all the antibiotics tested. We also reported the 
production of an Ambler class A P-lactamase GES-1 in a CC235 
isolate in Senegal. CC235 has been identified worldwide, in 
association with the carbapenemases VIM, IMP, FIM, and NDM 
[30]. This clonal complex falls into the CC /BURST Group 11 
(BGll) with a serotype Oil. GES- 1 -producing CC235 isolates 
have also been reported in Spain (Figure 3) [31]. Interestingly, all 
the CC235 isolates reported here [n = 8) were resistant to at least 
one antibiotic family. Finally, the last ESBL producer was the 
unrelated ST 1469 found in Senegal that carried WfloES-g 
(Figure 2). The gene WflcES-g has only been reported once in an 
isolate of P. aeruginosa from a patient hospitalized in France [32]. 
These data indicated that the presence of P-lactamases with an 
extended-spectrum in West and Central Africa is due to both the 
presence of successful international clonal complexes (CC233, 
CC235, and CC244) and that of distinct unrelated STs (ST1469). 
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Figure 3. Distribution of tKie STs and CCs of the Pseudomonas aeruginosa containing isolates producing either extended-spectrum p- 
lactamase or metallo-p-iactamase defined in the following references [5,28,45,46] on a dendrogram built with the data all known STs 
(n= 1595). ESBL- and MBL-producers were represented with red spots. See Figure 2 legend for the tree construction details. Pink zones surround STs 
that belong to the same clonal complex. 
doi:1 0.1 371 /journal.pone.01 07008.g003 



Resistance patterns of intercontinental clonal complexes 

Many epidemiological studies focused on multidrug-resistant 
isolates. Consequently, the vast majority of the isolates genotyped 
with MLST method are representatives of multi-drug resistant and 
outbreak STs. All microbiologists are encouraged to submit not 
only new allele sequences but also isolate information from 
already-known allele to the MLST databases so as it becomes a 
comprehensive resource for global epidemiology of this species. 
However, the submission of data for isolates belonging to already- 
known STs is not the rule. This constitutes one limitation of the 
databases. The spread of the multidrug-resistant intercontinental 
STs could be due either to their resistance pattern that favours 
tliem when the antibiotic pressure is high or to their intrinsic traits, 
independendy of the antibiotic resistance. A high-level resistance 
to major anti-pseudomonal agents can be the result of the 
acquisition of foreign resistance genes such as those encoding 
carbapenemase or of mutations that trigger chromosomal 
resistance mechanisms. To the best of our knowledge, the data 
published so far could hardly answer this question. One strong 
point of our study is the absence of selection biases of the clinical 
isolates. Hence, we genotyped every isolate of the collection 
whatever its resistance phenotype. All the CC233 and CC235 
isolates (w = 3 and 8, respectively) that we tested were multidrug- 
resistant, with the exception of two isolates resistant only to 
carboxypeniciUins. Cabot et al. confirmed the absence of these 
clonal complexes in a collection of multi-susceptible clinical 



isolates in Spain [33]. In contrast, all the isolates belonging to 
CC155 (w = 23), CC195 (n = 4), CC395 (n = 3) were fuUy 
susceptible to all the antipseudomonal compounds tested. How- 
ever, CC155 associated with 6/aviM-2 has been reported in Spain, 
and multi-drug resistant CC395 was countrywide-distributed in 
Hungary or epidemic in a French University Hospital [34,35]. 

CC244 isolates that carry 6/apER-i or 6/avEB-i have been 
reported in Poland and in Bulgaria [29,36]. We also found, in our 
collection, isolates producing the VIM-2 carbapenemase (w - 2) 
even though most of the isolates of this clonal complex (15 out of 
19) displayed a fuUy susceptible profile. Similarly, although the 
CC274 was also identified within multidrug-resistant isolates 
circulating in Spain [37], we found that the majority of the 
representatives of this clonal complex in our collection were of 
wild-type susceptibility pattern. 

Beside intercontinental clonal complexes, we observed the 
spread of fully susceptible STs in West and Central African 
countries. Hence, seven ST699 isolates were found in Ivory Coast, 
Nigeria, or Senegal and four ST 1464 isolates were found in 
Senegal. Altogether, these data could suggest that the spread of a 
ST among patients can be the result either (i) to its antibiotic 
resistance (e.g., CC235 and CC233) or (ii) to features independent 
from the resistance to antibiotics (e.g., CC155, CC244, CC195, 
CC274, and CC395). These latter features could favour human 
infection or colonisation directly with a better adaptation to the 
host (e.g. through a better adhesion to eukaryotic cells, higher 



PLOS ONE I www.plosone.org 



6 



September 2014 | Volume 9 | Issue 9 | el 07008 



Population Structure of Clinical P. aeruginosa from Africa 



North America 
94 isolates, 67 STs 



South America 
73 isolates, 17 STs 




Asia 

133 isolates, 109 STs 



Europe 

744 isolates, 497 STs 



••>' 



Oceania 

490 isolates, 225 STs 



\ 



Figure 4. Distribution of the STs of clinical isolates of Pseudomonas aeruginosa on a dendrogram built with the data all known STs 
(/7= 1595), broken down by their continent of origin. The geographical origin of each ST was extracted from the P. aeruginosa MLST database 
in July 2013 [11]. 

doi:1 0.1 371/journal.pone.01 07008.g004 



biofilm formation and virulence), or indirectly witli a higher 
frequency in environmental settings. Further studies are needed to 
clarify this point. 

Population structure of and acquisition of extended- 
spectrum p-lactamases by P. aeruginosa 

We sought to determine whether the STs of clinical Isolates of 
P. aeruginosa found in the West and Central Africa clustered on 
the global dendrogram of the 1595 currently defined STs 
(Figure 2). The analysis of MLST data relies on two strategies: 
allele- or nucleotide-based methods. The former is fast, simple, 
and very popular, but cannot distinguish between single-based 
changes in multiple loci versus multiple mutations in the same loci. 
It therefore disregards much of the evolutionary information 
contained at the nucleotide level. Nucleotide-based methods are 
more accurate and robust but require bioinformatics computer 
skills [38]. Figure 2 represents the ML tree produced by the 
alignment of the concatenated gene sequences of the currently 



known STs. Its bush-like structure confirmed that P. aeruginosa 
has a non-clonal structure [6-10]. The Figure 2 clearly shows the 
absence of any population substructures, confirming that P. 
aeruginosa evolves as a single cohesive genetic group. This, 
together with the large distance between the tree root and the 
representative of the close species P. fluorescens PfD-1 taken as 
outgroup confirmed that P. aeruginosa conformed to a clear, 
compact, and well defined species [39]. 

As the 80 STs identified our isolates were scattered throughout 
the ML tree, they were highly representative of the global P. 
aeruginosa population. This suggests no evidence of STs or clades 
specific to a region. We also spotted on the global ML tree the STs 
found on the five other inhabited continents (Figure 4) [11]. 
Although the number of STs submitted by continent is uneven 
(e.g., 497 in Europe and 17 in South America), it also suggests the 
poor association between a phylogroup and a continent. Using 
other methods (eBURST and calculation of genetic diversity 
index), Kidd et al. also showed no significant difference between 
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the genetic divc-rsity of STs found in Queensland (Australia) and 
the global P. aeruginosa ST dataset [10,11]. 

Standard phylogenetic methods assume a lack of recombination 
defined as the exchange of genetic information between two 
nucleotide sequences. Although the recombination play.s a major 
role in the evolutionary history of P. aeruginosa genome, the 
housekeeping genes used for MLST are presumably not subjected 
to these mutational events [40]. However, the phylogenetic tree 
based on these data should be analyzed with caution. In addition, 
MLST loci such as acs, guaA and mutL underlie a high mutation 
rate, particularly in isolates associated with chronic infection such 
as cystic fibrosis [41-43]. Consequently, it is likely that the genetic 
diversity of the species is lower than that calculated from the 
MLST data. 

Genetic diversity of P. aeruginosa producing |3- 
lactamases with an extended-spectrum 

Although Wombl or i/flESBL genes can be borne by many STs 
over the world, some successful international clonal complexes are 
frequently multidrug-resistant and their representatives are more 
prone to produce ESBLs and MBLs [5]. Members of these 
epidemic lineages have a serotype Oil and 012. The 012 and 
Oil lineages fall into the P. aeruginosa CC/BURST Group 4 
(BG4) and BGl 1 and embrace STs that cluster tightly in the ML 
tree (Figure 3). We confirm here that MBL and ESBL-encoding 
genes are mostly borne by successful international CC233, 
CC235, and CC244 as previously observed in other countries 
[28,30,44]. The graphical representation on the global ML tree of 
all the STs defined so far as potential MBL- or ESBL-producers 
suggested that the acquisition of this resistance enzyme has a very 
limited association with a genotype [5,28,45,46] (Figure 3). 

In conclusion, we found that the resistance rates of clinical 
isolates to antibiotics remained low in Senegal and Central African 
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